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239  Fronczak  Hall 

State  University  of  New  York  at  Buffalo 
Buffalo,  New  York  14260 


Abstract 

In  order  to  calculate  different  electronic  states  in  polyatomic  rare 
gas-halogen  systems,  including  ionic  states  with  charge  transfer,  a 
semiempirical  diatomics- in- ionic-systems  (DIIS)  method  is  presented.  In 
this  method  the  Hamiltonian  matrix  elements  are  expressed  in  terms  of 
diatomic  potentials,  except  for  those  which  are  responsible  for  coupling 
between  the  ionic  and  neutral  states.  The  coupling  matrix  elements  are 
determined  by  fitting  the  diatomic  transition  moments.  Calculations  are 
performed  for  Xe2Cl,  Xe^Cl  and  Xe^Cl.  The  results  reveal  at  least  three 
quasistable  ionic  Xe^Cl  molecules.  Such  molecules  can  attract  more  Xe 
atoms,  which  are  binded  mainly  by  polarization  forces. 
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I .  Introduction 

Interest  in  rare  gas-halogen  systems  has  been  increasing  due  to  their 

1-12 

spectroscopic,  and  in  particular  lasing,  properties.  The  lasing 

properties  are  associated  with  the  formation  of  quasistable  ionic  molecules 
R+X  and  R*X  ,  where  R  and  X  are  the  rare  gas  and  halogen  atoms, 
respectively.  These  ionic  molecules  are  obtained  in  rare  gas-halogen  gas 

mixtures  either  by  electron  impact*  or  by  optical,  most  often  laser, 

5  9  11  12  +  - 

irradiation.  *  *  *  Recently  the  ionic  rare  gas-halogen  molecules  Xe^Cl 

13  14 

were  detected  also  in  rare  gas  solids  doped  by  Cl^  and  HC1.  *  Rare  gas- 

halogen  systems  are  of  interest  not  only  for  the  configurations  where  the 

quasistable  ionic  R+X  and  R*X  molecules  are  formed,  but  also  for 

15-19 

configurations  of  van  der  Waals  RX,  RX2  or  RXH  complexes. 

The  theoretical  consideration  of  the  rare  gas-halogen  systems  involves 

a  knowledge  of  potential  curves  and  potential  energy  surfaces  (PES).  The 

calculation  of  these  PES  is  the  aim  of  the  present  paper.  For  the  sake  of 

generalization,  we  shall  consider  here  arbitrary  RnX  systems  where  n  may  be 

larger  than  2  and  consequently  the  number  of  atoms  larger  than  3.  Although 

the  RX  and  R2X  systems  have  been  mostly  the  subject  of  study,  the  R^X 

systems  with  n  >  2  can  also  be  of  interest  due  to  the  possibility  of  the 

formation  of  cluster-like  molecules  or  large  van  der  Waals  complexes.  The 

R*X  (n  >  2)  excited  systems  can  also  be  formed  in  rare  gas  solids  doped  by 
20 

halogen  molecules. 

In  the  case  of  diatomic  RX  (n  *  1)  systems  some  of  the  potential  curves 

for  different  electronic  states  are  well  known  from  ab  initio 
21-23 

calculations.  The  RX  potential  curves,  mostly  for  the  ground  state, 

24 

were  studied  also  by  applying  model  approach  and  by  using  experimental 
data.^’*^’^  The  information  concerning  the  triatomic  R2X  systems,  not  to 


mention  the  R^X  systems  for  n  >  2,  is  much  more  scarce.  The  experimental 

data  of  the  R£X  systems  are  limited  mainly  to  the  transition 
7  9  11  12 

energies.  *  *  *  The  all-electron  ab  Initio  calculations  were  performed 

27 

for  the  Ar^F  and  K^F  systems  only.  The  ab  initio  calculation  of  the 

Xe2Cl  system  which  contains  many  more  electrons  was  performed  by  using  a 
28 

simplified  model.  The  R^X  systems  were  also  calculated  semiempirically  by 

29  30 

using  the  diatomics- in-molecules  (DIM)  method.  * 

The  advantages  of  the  semiempirical  approach  in  the  PES  calculations  of 

polyatomic  systems  containing  a  large  number  of  electrons  are  well  known. 

37 

In  particular,  the  DIM  method  is  often  used  in  the  PES  calculations. 

Unfortunately,  in  the  case  where  ionic  states  are  involved,  the  application 

of  the  DIM  method  becomes  much  more  complicated  if  the  coupling  between  the 

neutral  and  ionic  states  is  taken  into  account.  The  difficulties  arise 

largely  due  to  the  "zero  overlap  of  atomic  orbitals"  approximation  which  is 

usually  used  in  the  DIM  method.  According  to  this  approximation,  the 

overlap  between  atomic  orbitals  belonging  to  different  atoms  is  neglected. 

Consequently  the  off-diagonal  matrix  elements  between  any  electronic 

configurations  with  different  numbers  of  electrons  on  the  same  atom,  like 

the  neutral  and  ionic  states,  are  equal  to  zero,  such  that  the  states  are 

uncoupled.  This  shortcoming  of  the  usual  DIM  method  can  be  overcome  by 

32-34 

introducing  additional  assumptions  and  new  parameters.  However,  such 

an  approach  increases  the  complexity  of  the  calculation  and  makes  the 

29  30 

results  more  arbitrary.  In  DIM  calculations  of  the  R^X  systems,  *  the 
usual  version  of  the  method  was  used  with  the  off-diagonal  matrix  elements 
between  neutral  and  ionic  states  equal  to  zero  and  the  independent 
consideration  of  the  R2X  and  R^X  states. 


The  influence  of  the  coupling  between  the  neutral  and  ionic  states  of 

the  R^X  systems  is  not  quite  clear.  This  problem  was  studied  only  in  the 

22 

connection  with  ab  initio  calculations  of  the  XeX  diatomic  systems. 

According  to  this  study,  the  coupling  is  of  the  order  of  10%  for  the 

equilibrium  XeX  van  der  Waals  distance  and  decreases  slowly  as  the 

interatomic  distance  increases.  It  is  difficult  to  estimate  the  influence 

of  the  coupling  on  different  properties  of  the  RX  systems,  although  it 

obviously  reduces  size  of  the  van  der  Waals  RX  molecule.  For  example,  the 

35 

van  der  Waals  radii  of  the  Xe  and  Cl  atoms  are  2.2  A  and  1.8  A, 

respectively,  and  while  the  XeCl  van  der  Waals  distance  is  expected  to  be 

4.0  A,  in  reality  it  is  only  3.23  A.  Because  the  coupling  may  be 

different  for  different  rare  gas-halogen  systems,  the  use  of  diatomic  RX 

potential  curves  in  the  semiempirical  calculations  of  polyatomic  systems 

(DIM,  for  example)  may  lead  to  confusing  results.  Thus,  the  coupling  and 

consequently  the  interaction  are  expected  to  be  small  in  the  case  of  a  RX£ 

van  der  Waals  complex  due,  partly,  to  lower  electron  affinity  of  the  X2 

20 

molecule  compared  to  the  X  atom,  which  is  confirmed  by  the  spectroscopic 

36 

properties  of  X2  molecules  trapped  in  a  rare  gas  matrix.  However, 
neglecting  the  change  in  the  coupling  and  using  the  empirical  RX  interaction 
in  the  DIM  calculation  of  the  van  der  Waals  RX2  molecule,  one  would  obtain 
the  RX2  interaction  as  strong,  roughly,  as  in  the  R-X  case,  in  contradiction 
with  the  experiment.  The  coupling  between  neutral  and  ionic  states  must  be 
taken  into  account  if  one  wants  to  calculate  the  off-diagonal  (transition) 
dipole  moments. 

We  shall  present  in  this  paper  a  semiempirical  method  which  is  designed 
for  a  calculation  of  rare  gas-halogen  systems  RnX  consisting  of  an  arbitrary 
number  of  rare  gas  atoms  and  one  halogen  atom.  Introducing  some 


3 


modifications,  the  method  can  be  extended  to  R^HX  systems.  The  method  takes 


into  account  in  a  straightforward  way  the  coupling  between  the  neutral  R^X 


and  the  ionic  RfiX  configurations,  as  well  as  the  positive  charge 


delocalization  in  the  ionic  R  X  states.  The  method  makes  use  of  diatomic 

n 


potentials  as  initial  data  for  a  calculation,  but  in  a  different  way  than 


the  DIM  method.  A  main  feature  of  the  method  is  the  consideration  of  the 


systems  with  ionic  states,  and  hence  we  call  it  the  "diatomics- in- ionic- 


systems"  (DIIS)  method. 


The  spin-orbit  coupling  in  the  Cl  atom  and  Xe  ion  will  be  not  taken 


into  account  directly.  The  shift  of  energy  levels  resulting  from  this 


coupling  can  always  be  incorporated  indirectly  by  using  proper  empirical  or 


ab  initio  potentials.  The  main  disadvantage  of  the  neglect  of  spin-orbit 


coupling  is  the  loss  of  some  of  the  excited  states.  The  direct  inclusion  of 


spin-orbit  coupling  does  not  present  theoretical  obstacles.  However,  it 


would  make  the  Hamiltonian  matrix  elements  much  more  complicated  and  would 


increase  significantly  the  order  of  the  matrix.  The  latter  is  of  primary 


concern  for  large  systems,  like  rare  gas  solids  with  halogens. 


The  theoretical  consideration  of  the  DIIS  method  is  given  in  Section 


II.  Section  III  presents  the  equations  for  the  diatomic  terms  which 


determine  the  matrix  elements  of  the  method.  The  diatomic  terms  for  the 


XenCl  systems  are  given  in  Section  IV.  Section  V  deals  with  the  calculation 


of  the  Xe^Cl  systems,  namely  Xe^Cl,  Xe^Cl  and  Xe^Cl*  and  Section  VI  is  the 


Conclusion. 


II.  Diatomics-in-ionic-svstems  method 


A  polyatomic  rare  gas-halogen  system  RfiX  will  be  presented  here  as 


"A-  Rj _^Xj  (J  =  n  +  1),  where  the  R^  atoms  have  a  closed  shell  in  the 


neutral  state  and  P  symmetry  shell  in  the  ionized  state  R,,  whereas  the  X 


•'»  -I.-1*  v*stU  .VoTv 


M2 


atom  has  a  P  symmetry  shell  in  the  neutral  state  and  a  closed  shell  when 

becoming  the  negative  ion  XT.  The  R.  atoms  need  not  be  just  the  rare  gas 

atoms  but  also  H  atoms,  with  the  difference  that  in  the  ionized  state  H+ 

they  are  also  spherically  symmetric.  The  systems  with  H  atoms  can  be 

treated  as  a  particular  case  which  we  shall  not  consider  separately, 

suggesting  that  all  R*  ions  have  a  P  symmetry  shell. 

The  wave  function  of  the  polyatomic  system  in  a  fixed  electronic 

configuration  will  be  presented  here  as  an  antisymmetrized  product  of  atomic 
37 

group  functions,  which  are  assumed  to  be  exact  many-electronic  wave 
functions  of  individual  atoms. 


*  =  A  n  X.  ,  (1) 

1  J  i  l 

where  A  is  the  operator  which  antisymmetrized  electrons  belonging  to 
different  group  functions  xim  and  are  the  atomic  state  indexes.  The 

atooic  group  functions  in  the  polyatomic  wave  function  (1)  are  not 

orthogonal  one  to  another  which  makes  the  energy  expression  very 
38,39 


complicated. 


The  consideration  is  significantly  simplified  by  using 


the  zero  overlap  of  atomic  orbitals  (ZOAO)  approximation  which  is  very 

31 

comaon  in  semieropirical  methods,  in  particular  in  the  DIM  method.  The 

ZOAO  approximation  allows  one  to  omit  the  antisymmetrization  operator  A, 

thus  rendering  the  polyatomic  wave  function  as  a  simple  product  of  atomic 

group  functions.  The  mistake  arising  from  the  neglect  of  the  overlap  is 

39 

proportional  to  the  square  of  the  interatomic  overlap  integrals  and  has 

40 

been  shown  to  be  small,  at  least  in  DIM  calculations. 

In  the  accepted  approximations,  the  wave  function  of  a  neutral 
Rj_jXj  electronic  configuration  is  equal  to  the  product  of  S 


5 


symoetrical  (closed  shell)  group  functions  of  the 

symmetrical  group  function  of  the  XT  (halogen)  atom, 

jm  j 


atoms  and  the  P 


Jm 


m 


1,2,3 


(2) 


where  the  index  m  indicates  the  orientation  of  the  Xj^  group  function  along 
the  x,y,z-axes.  In  an  ionic  state,  one  of  the  J-l  rare  gas  atoms  (R*)  has  a 
positive  charge  and  is  distinguished  by  its  X^j^  group  function  orientation, 
whereas  the  Xj  ^  group  function  of  the  halogen  ion  Xj  has  S  symmetry, 

J-l 


4>.  =  n  x-^^X^+^xl  \  i=l»2...J-l;  m=l,2,3;  j*i. 

lm  .  Aj  Aim  AJ  J 

J 


(3) 


The  upper  indexes  (0),  (+)  and  (-)  indicate  in  Eqs.  ( 2 ) - ( 3 )  and  below  the 
charge  of  an  atom.  It  is  important  to  emphasize  that  Eqs.  ( 1 ) - ( 3 )  do  not 
necessarily  imply  the  one-electron  approximation,  since  the  group  functions 
are  arbitrary  and  can  include  correlation  between  electrons  belonging  to  the 
same  atom.  Using  the  group  functions  product,  we  neglect  only  the 
interatomic  correlation  which  describes,  in  particular,  the  van  der  Waals 
attraction.  However,  in  the  framework  of  the  semiempirical  approach,  the 
van  der  Waals  interaction  can  be  taken  into  account  indirectly  by  using 
empirical  diatomic  potentials. 

The  three  functions  (Eq.(2)]  of  the  neutral  configurations  and  the 
3(J-1)  functions  (Eq.(3)]  of  the  ionic  configurations  form  the  polyatomic 
basis  set.  The  linear  combination  of  these  functions  forms  the  wave 
function  of  the  system 


* 


J  3 


l  l 

1=1  m*l 


(4) 
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m 


>1 


The  energy  levels  of  this  wave  function  are  obtained  now  in  the  usual  way  by 


forming  the  Hamiltonian  matrix  of  order  3J*3J,  where  J  is  the  number  of 


atoms  in  the  system. 


Due  to  the  lack  of  interatomic  electron  permutations  in  the  polyatomic 


functions  of  Eqs.(2)  and  (3),  the  Hamiltonian  of  the  system  can  be 


partitioned  into  atomic  and  interatomic  terms.  This  partitioning  depends  on 


the  representation  of  the  polyatomic  wave  function  and  changes  from  one 


matrix  element  to  another.  The  Hamiltonian  which  is  used  for  the  evolution 


of  the  matrix  elements  for  the  neutral  states  has  the  form 


J-2  J-l 


■[  W*!  l  hW VV +  W + 1  vvv- 

j  Ji  j 


where  the  ft.  and  h  terms  with  one  index  include  all  energy  operators  of  the 
J  ^ 


electrons  belonging  to  one  atom,  the  h.  .  and  h.T  terms  with  two  indexes 

J1J2 


include  the  interactions  between  electrons  and  nuclei  belonging  to  different 


atoms,  and  indicates  the  number  of  electrons  in  atom  j.  We  will  present 


the  energy  of  the  system  relative  to  the  energy  of  separated  on  infinity 


atoms 


*.  *  l  EJ0)  +  E$0)> 


where  E*^  and  are  the  energies  of  the  rare  gas  and  halogen  atoms, 

J  J 


respectively.  Using  the  expressions  (2),  (5)  and  (6),  one  obtains  the 


following  expression  for  the  diagonal  matrix  elements  of  the  neutral  states: 


Jm ,  Jm 


fdv*T  ft*. 

J  Jm  Jm 


J-2  J-l 


l  l  u<vj2  + 1 

ji  J2(>Ji)  j 


In  Eq.  (7),  U''?'".  and  IT  are  the  potentials  of  the  interaction  between 

Jl*J2 

two  rare  gas  atoms  and  between  a  rare  gas  atom  and  a  halogen  atom, 
respectively: 


u<°°>,  -  h.  .  [x<°>r 

Jl » J2  J  J  j  ^1^2  J2 

■  Ko)i*  v-£>l2- 


In  the  expression  for  the  off-diagonal  matrix  elements  of  the  neutral 

states,  all  integrals  with  one-atom  operators  and  and  with  the 

interatomic  terms  h.  .  vanish  either  due  to  the  difference  in  the  symmetry 

J1J2 

of  functions  with  different  m  or  due  to  the  ZOAO  approximation.  Only  the 

h  terms  give  a  non-zero  contribution, 

J  J 


Ht  t  =  fdv$  ft*  =  V  fdvlx^ 
Jmj.Jn^  J  Jm1  Jm2  £  J 


(»a 

j  1  jJxJn  jJm,' 


Let  us  consider  now  the  diagonal  matrix  elements  for  the  ionic  states 

[ Eq - ( 3 ) ] .  The  Hamiltonian  for  the  function  which  describes  the  state 
^  lm 

with  a  positive  charge  on  the  i-th  rare  gas  atom  and  a  negative  charge  on 
the  J-th  (halogen)  atom  has  the  form 

J-l  J-2  J-l  J-l 

* ■  I  W + 1  l  V<VV +  + 1  Vv'+v* 

J  j2<>J1J  1212  j 

J-l 

fijCvj+l)  +  £  hjjCVj+Vj+l)  +  hiJ^vi+vj^»  *  i’ 

j 

This  form  includes  terms  which  are  responsible  for  the  polarization  of 
neutral  atoms  by  the  ions  R*  and  Xj.  It  is  more  convenient  to  determine  the 
energy  of  polarization  in  the  electrostatic  approximation,  separating  it  as 


*  V  i 


'M'S  Vi 


t't- 

iv 


a 


I 


*y.i 


W 

>•* 

I 

'§ 


*  *.*1 


t’il 


a  special  term 

H.  .  =  8.  .  +  £, • 

im,im  im,im  i 

Considering  atoms  and  ions  as  non-overlaping  spheres,  one  obtains 


1. 


p  A,  ,  rv  - 

-  -V  C ,(-U-  -  — 
1  L  j  o  ^  R 


R..  R,. 

1J  Jj 


3>  +  a.,  j*i, 


C.  =eV/2, 

where  a.  are  the  rare  gas  atom  polarizabilities,  and  are  the 

vectors  between  the  corresponding  atoms,  and  includes  the  dipole  dipole 

interaction.  The  matrix  elements  are  obtained  by  using  the 

_  ,_4 

Hamiltonian  without  polarization  terms  -C^/R^: 
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Using  the  expressions  (3),  (6)  and  (101),  one  obtains 
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The  off-diagonal  matrix  elements  between  the  ionic  states  with  the  same 
localization  of  the  positive  charge  but  different  orientation  m  are  obtained 


in  the  same  way  as  the  matrix  elements  in  Eq.(9), 
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The  off-diagonal  matrix  elements  between  the  ionic  states  with 
different  localization  of  the  positive  charge  and  between  the  neutral  and 
ionic  states  are  also  obtained  in  a  similar  way, 
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where  H.  .  (v.  +v.  -1)  is  the  Hamiltonian  of  the  charged  (R.  R.  )+  diatomic 
1112  11  x2  *1  X2 

fragment  and  fi  .(v  +v  )  is  the  Hamiltonian  of  the  R.X-  diatomic  fragment. 

The  exchange  integral  (18)  presents  the  valence  interaction  between  two  rare 
gas  atoms  with  a  common  +1  charge.  The  exchange  integral  (19)  is 
responsible  for  the  coupling  between  the  neutral  and  ionic  states  of  the 
rare  gas-halogen  diatomic  systems. 

The  matrix  elements  obtained  above  are  expressed  by  diatomic  terms 
which  are  determined  by  atomic  group  functions  oriented,  when  of  P  symmetry, 
along  Cartesian  (x,y,z)  axes.  In  order  to  express  the  matrix  elements  by 
the  diatomic  terms  of  usual  E  and  II  orientations,  one  needs  to  transform  the 
P  symmetrical  group  functions  to  a  new  set  of  functions 


xim  \  ^ij,mk^ik,j' 


(20) 


k-1 


where  the  indices  mal,2,3  indicate  the  x-,  y-  and  z-orientations, 

respectively,  <|»..  is  the  E  group  function  of  the  i-th  atom  oriented  along 
ii »  j 


the  1?.  .  vector  connecting  the  two  atoms,  and  ij/  _  and  <|*  -  are  the  group 
ij  » j  » J 

functions  of  the  i-th  atom  oriented  perpendicular  to  the  -vector.  We 
will  choose  one  of  these  two  perpendicular  orientations  (n)  as  located  in 
the  plane  formed  by  the  -vector  and  the  x-axis.  The  second  orientation 
which  is  perpendicular  to  this  plane  will  be  denoted  by  S.  The  3*3  mlf- 


matrix  has  the  form 


S.,  , ,  =  S.  .  _  =  (x.-x.)/R..,  S 
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-(xj-xi)(zj-zi)/(Rijr.j),  S.j33  =  S.j>zfi  «  -(yj*yi)/r.j,  (21") 

where  R^  is  the  interatomic  distance  and  r^  is  the  projection  of  the  R^ 
vector  on  the  yz-plane. 
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Using  the  transformation  (21)  we  will  express  the  matrix  elements  (7),  (9), 
(15),  (17),  (18)  and  (19)  by  Z  and  n  diatomic  terms, 
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The  diatomic  terms  V.  .  and  V..  in  the  off-diagonal  matrix  elements  (27) 

V2 


and  (28)  are  expressed  by  the  integrals  (18)  and  (19),  respectively,  with 


the  only  difference  that  instead  of  x,y,z-oriented  functions  x>  the  Z,H- 


oriented  functions  i|»  of  Eq.  (20)  are  used.  The  diatomic  U  terms  in  the 


matrix  elements  ( 23 ) - ( 26 )  will  be  called  Coulombic  or  diabatic  potentials  as 


they  do  not  include  exchange  interactions  presented  by  the  off-diagonal 


Y*  V 


matrix  elements  (27)-(28).  We  have  to  note  that  the  Coulombic  terms  U  can 


include  in  a  semiempirical  way  the  van  der  Waals  interaction  which  is  not  of 
Coulombic  character. 


The  3Jx3J  matrix  H  provides  3J  eigenvalues  and  3J  eigenvectors  c 

.00 


(k) 

im 


of  the  wave  function  of  Eq.(4).  Using  the  eigenvectors  c.  ,  it  is  possible 

im 

to  calculate  both  static  and  transition  moments.  In  the  ZOAO  approximation 

the  static  dipole  moment  of  any  state  is  completely  determined  by  the 

contribution  of  the  ionic  configuration.  The  x-component  of  the  dipole 

moment  u,  is 
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where  e  is  the  electron  charge  and  and  are  the  radius  vectors  of  the 

J  X 

ions.  The  transition  moment  between  a  state  with  largely  ionic 
configuration  and  a  state  with  largely  neutral  configuration  is 
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III.  Diatomic  Coulombic  and  exchange  terms 

In  the  previous  section  the  matrix  elements  of  a  polyatomic  system  were 
expressed  by  diatomic  Coulombic  (U)  and  exchange  (V)  terms.  We  shall  now 
express  these  Coulombic  and  exchange  terms  by  diatomic  potentials  which  are 
assumed  to  be  known.  For  the  sake  of  simplicity,  our  study  will  be 


restricted  to  systems  with  identical  rare  gas  atoms  R^.  There  are  four 
kinds  of  diatomic  fragments:  R^,  R^  *  RX  and  RX.  The  expressions  for  the 
energy  of  any  diatomic  fragment  will  be  obtained  from  the  matrix  elements  of 
the  previous  section  by  removing  all  atoms  to  infinity  except  those  which 
belong  to  this  fragment. 

A.  R-R  interaction 

In  the  case  of  the  R^R^  diatomic  fragment,  the  ground  state  energy  is 

equal  to  the  diagonal  matrix  element  [ Eq . ( 23 ) ]  with  one  non-zero  term  only, 

namely  It  follows  that  the  diatomic  term  is  identical  to  the 

J  ij  ij 

R2  ground-state  I  potential 

U?00)  =  I.  (33) 

ij 

B.  (R-R)+  interaction 

The  energy  of  the  diatomic  fragment  (R^R^)+  is  determined  by  a  2*2 
matrix  with  non-zero  off-diagonal  matrix  elements.  In  the  case  of  the  I- 
oriented  R+  wave  function,  the  diagonal  matrix  element  (24)  and  the  off- 
diagonal  matrix  element  (27)  become 

Hn  ’  "22  *  V  V  *  cXj  *  V  V  <34) 

H12  ‘  <34,) 


where  -C^/R^j  is  the  polarization  energy  e.  The  matrix  elements  (34)  and 
(34')  provide  a  ground-state  potential  and  an  excited  state  potential 
Considering  these  potentials  relative  to  the  energy  of  separated  ion  R+  and 


neutral  atom  R,  one  obtains 

uij,s  1(IU  +  V*  °ij  ,z  “ij.i  +  VR1J 

''ij.!  -  »<£u  -  V' 


Similar  expressions  are  valid  also  for  the  n  states,  with  the  only 


(37) 


D.  R-X  interaction 

Because  of  the  ionic-neutral  coupling,  the  R-X  diatomic  potentials  are 
not  connected  explicitly  with  the  Coulomb  terms  U^j,  like  in  the  previous 
cases,  but  rather  are  obtained  as  the  eigenvalues  of  a  2*2  matrix  with  the 
diagonal  elements  expressed  by  the  Coulombic  terms  Ih  j  and  the  off-diagonal 
elements  equal  to  the  exchange  term  V^j.  In  the  case  of  the  E-orientation, 
the  matrix  elements,  derived  from  Eqs.  (23), (24)  and  (28),  are 
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The  matrix  (40)  allows  us  to  express  the  Coulombic  terms  an(j  ^  in 

terms  of  its  eigenvalues,  i.e.,  the  RX  ground-state  (1Z)  and  excited-state 

(2E)  potentials  which  are  assumed  to  be  known: 
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3  -  V(2E-1E)2  -  4V2  .  (43) 

Similar  expressions  are  valid  for  the  n-oriented  potentials  ujj^^  and  UR+  \ 

Equations  (4l)-(42)  exclude  from  the  R-X  potentials  the  ionic-neutral 

copupling  term  V^j,  which  must  be  found  independently.  Before  estimating 

this  term,  we  shall  consider  the  influence  of  the  ionic-neutral  coupling  on 

the  R-X  properties.  Due  to  this  coupling,  the  ground  state  (IE)  of  the  rare 

gas-halogen  molecule  is  a  mix  of  "pure"  neutral  and  ionic  states  whose 

energies  are  determined  by  the  diagonal  matrix  elements  (40)  and  (40'). 

Because  the  ionization  potential  of  the  rare  gas  atom  IR  is  much  higher  than 

the  electron  affinity  of  the  halogen  atom  A^,  the  ^22 ~^11  difference  *s 

large.  It  results  in  relatively  weak  coupling,  at  least  for  distances 

20  22 

exceeding  the  van  der  Waals  equilibrium  distance.  ’  However,  this 
coupling  significantly  affects  some  of  the  properties  of  the  rare  gas- 
halogen  systems.  These  properties  can  be  used  to  find  the  coupling  terms  by 
semiempirical  fitting.  In  the  present  calculation  the  coupling  term  is 
found  by  fitting  the  RX-to-R+X  transition  moment  which  depends  strongly  on 
V^j  and  vanishes  when  is  zero. 

IV.  Diatomic  terms  of  the  Xe.  Cl  systems 

For  the  XenCl  systems  the  diatomic  fragments  are:  Xe-Xe,  (Xe-Xe)+, 
Xe-Cl  and  Xe-Cl.  The  Coulombic  potentials  and  the  exchange  terms  of  all  of 
these  diatomic  fragments,  except  Xe-Cl,  are  expressed  by  diatomic  E  or  II 
potentials  (see  Eqs.  (33)-(39)].  In  the  case  of  the  Xe-Cl  diatomic 
fragment,  as  it  was  shown  in  the  previous  section,  the  Coulombic  potentials 
U  [Eqs. (41 ) - ( 42 ) ]  are  expressed  not  only  by  diatomic  potentials  but  also  by 
the  exchange  (coupling)  term  V,  which  has  to  be  found  independently. 


The  diatomic  potentials  can  be  taken  either  from  experimental  studies 
or  from  ab  initio  calculations.  In  the  present  calculation,  when  possible, 
empirical  potentials  were  used.  The  diatomic  potentials  were  expressed 
analytically  by  the  equation 

U  -  Uq  +  Aexp( -aR]  -  C6/R*  -  CjRU  -  C^R,  (44) 

where  R  is  the  interatomic  distance,  Uq  is  the  energy  of  two  separated 
atoms  or  ions,  is  the  van  der  Waals  coefficient,  is  the  polarization 
coefficient  (Eq.(13)],  and  Cj  *  e  is  the  coefficient  of  the  electrostatic 
attraction  between  two  opposite  charges.  In  the  case  of  a  pair  of  neutral 
atoms,  only  the  repulsive  and  van  der  Waals  terms  remain  since  *  0: 

U  *  U  +  Aexp(-oR]  -  C, /R6  (44') 

o  o 

The  0  potential  ( Eq . ( 16 ) J  of  an  ion-atom  pair  is  also  expressed  by  Eq. 

(44'),  since  *  0  and  the  polarization  term  is  excluded.  For  a  better 
fitting,  some  of  the  potentials  can  include  two  exponential  functions  or  the 

g 

-Cg/R  term.  The  exchange  terms  V  are  expressed  by  an  exponential  function 

V  *  Aexp(-aR).  (45) 

The  parameters  of  Eqs.  (44)  and  (45)  are  presented  in  Table  1.  The  sources 
of  the  potentials  are  given  below. 

A.  Xe-Xe  potential 

The  parameters  of  the  expression  (44')  for  the  ground  state  Z 
potential  were  obtained  by  the  interpolation  of  the  experimental  data  of  the 
Xe-Xe  interaction.4^ 

B.  (Xe-Xe )*  potentials  and  exchange  terms 

The  ab  initio  Xet  potential  curves  Z  ,  I  ,  II  ,  II  without  spin-orbit 
- 2  y  g  u'  u’  g  Y 
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coupling  were  used  to  calculate  the  Coulombic  potentials  U  (Eqs. (35)  and 


(37)]  and  the  exchange  terms  V  (Eqs.  (36)  and  (38)].  These  potentials  are 


preferable  as  they  provide  the  ground-state  potential  well  (D  =  1.08eV,  R 

e  e 

*  3.22  A)  which  is  more  reliable  from  the  point  of  view  of  experimental 
A  A 

results  than  that  obtained  in  Ref.  43  by  including  the  spin-orbit 
coupling.  The  (Xe-Xe)+  diatomic  terms,  U  and  V,  were  interpolated  by  the 
analytical  expressions,  Eqs.(44)  and  (45). 

C.  Xe-Cl  potential 

The  Xe-Ci  Z  potential  was  obtained  by  interpolation  of  experimental 
45 

data.  In  order  to  obtain  better  fitting  for  short  distances,  the 
repulsive  term  is  described  by  two  exponential  functions. 

D.  Xe-Cl  potentials  and  coupling  terms  (Fig.  1) 

In  order  to  find  the  Xe-Cl  Coulombic  terms  and  \  one  needs 

the  Xe-Cl  potentials  and  the  coupling  term  V  [see  the  Eqs.  (4l)-(42)].  The 

potential  curves  of  the  ground  state  II  (X)  and  the  excited  ionic  state  2Z 

26 

(B)  were  obtained  by  the  interpolation  of  experimental  data.  The 

potential  curve  of  the  excited  ionic  211  (C)  state  is  similar  to  the  2Z  (B) 
22 

potential  curve.  The  equilibrium  point  energy  of  2n  was  found 

46 

experimentally  to  be  almost  the  same  as  in  the  2Z  state  or  to  be  lower  by 
8  47 

only  0.016  eV  or  0.028  eV.  The  211  equilibrium  distance  was  found  to  be 
.47  * 

Re  ■  3,074  A,  which  is  larger  by  -0.06  A  than  that  in  the  2Z  state.  The 
211  potential  curve  used  in  this  paper  has  a  minimum  at  the  point  Re  *  3.074 
A  and  minimum  energy  value  lower  than  the  2Z  state  by  0.02  eV,  which  is  a 
compromise  of  different  experimental  results. -j^e  p0iari2ation 


coefficient  C^  for  the  211  potential  is  assumed  to  be  the  same  as  for  the  2Z 
potential  curve. 

The  potential  curve  of  the  first  excited  in  (A)  state  was  found  to  have 
a  minimum  at  the  point  R^  ■  4.1  A.  The  in  potential  curve  used  in  the 

present  paper  fits  this  minimum  point  and  the  energy  of  transition  from  the 


equilibrium  point  of  the  211  state  to  the  repulsive  III  state,  which  is  equal 
to  3.59  eV  (345  nm).49 

By  using  empirical  Xe-Cl  potentials  we  take  into  account  indirectly  the 

spin-orbit  coupling,  but  only  as  a  shift  of  the  diatomic  terms  in  the 

diagonal  matrix  elements.  Any  influence  of  the  interaction  of  the  Xe-Cl 

fragment  with  other  atoms  on  the  spin-orbit  coupling  is  neglected.  Besides, 

the  consideration  is  restricted  to  states  which  correspond  asymptotically  to 

the  Xe+  and  Cl  atoms  in  their  state* 

As  mentioned  in  the  previous  section,  the  coupling  terms  V^.  and 

which  are  responsible  for  the  mixing  of  the  neutral  and  ionic  states  can  be 

considered  as  semiempirical  parameters.  In  the  present  calculation  they 

were  determined  by  fitting  the  XeCl  transition  moments.  The  ab  initio 

22 

transition  moments  were  found  as  functions  of  interatomic  distances.  The 

empirical  transition  moments  can  be  found  at  the  equilibrium  points  of  the 

22 

ionic  states  by  using  the  radiative  lifetime  , 

77-9<V  V'3/2(W1/2-  Ek  >  Et  (46) 

where  y  is  expressed  in  Debyes  (D),  E  in  eV,  and  t  in  nsec.  The  lifetimes 

of  the  B->X  and  C-+A  (Fig.  1)  transitions  are  equal  to  11.1  ns  and  131  ns, 

respectively.  Substituting  these  values  in  the  expression  (46),  as  well 

as  the  energy  differences  (4.02  eV  for  B-+X  and  3.59  eV  for  C-*A  transition), 

one  obtains  y^.  *  2.9  D  for  the  B->X  transition  moment  and  y^  *  1.0  D  for  the 

C-*A  transition  moment.  These  values  were  fitted  by  the  coupling  terms  V^. 

and  Vjj  presented  as  exponential  functions  [Eq.(45)J.  The  coefficients  a  of 

these  exponential  functions  were  determined  by  using  the  ab  initio 

22 

transition  moments.  The  parameters  of  the  coupling  functions  are  given  in 
Table  1.  The  dependence  of  the  B-X  and  C-A  transition  moments  on  the  Xe-Cl 
distance  is  demonstrated  in  Fig.  2a. 


Table  1.  The  dependence  of  the  B-X  and  C-A  transition  moments  on  the  Xe-Cl 
distance  is  demonstrated  in  Fig.  2a. 


The  effect  of  the  coupling  on  the  Xe-Cl  interaction  is  demonstrated  in 
Fig.  2b  by  presenting  the  coupling  energy,  which  is  equal  to  the  difference 
between  the  Coulorabic  potential  and  the  ground-state  potential  IE,  and 

the  contribution  (in  X)  of  the  ionic  Xe+Cl  configuration  to  the  ground 
state.  Both  the  coupling  energy  and  the  ionic  configuration  contribution 
decrease  exponentially  as  the  internuclear  distance  increases.  At  the  point 
of  the  van  der  Waals  equilibrium  (3.23  A),  the  coupling  energy  is  0.10  eV, 
almost  equal  to  the  energy  of  the  van  der  Waals  attraction  (0.12  eV).  The 
contribution  of  the  ionic  configuration  at  this  point  is  2.4%. 

V.  Results  of  calculations 

The  semiempirical  DIIS  calculations  were  performed  for  the  following 

Xe^Cl  systems:  Xe2Cl,  Xe3Cl  and  Xe^Cl.  The  Xe2Cl  system  has  been  studied 
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both  experimentally  '  *  ’  and  theoretically.  *  In  the  present 

30 

calculation,  which  unlike  another  semiempirical  calculation  takes  into 
account  directly  the  neutral-ionic  coupling,  the  transition  moments  were 
obtained  as  well.  The  Xe^Cl  and  Xe^Cl  systems  are  considered  here  in  order 
to  check  the  possibility  of  the  formation  of  more  stable  ionic  molecules 
than  the  Xe2Cl  molecule. 

A.  Xe^Cl  system 

According  to  the  results  of  the  present  calculation,  the  quasistable 

4-  -  2 

Xe2Cl  molecule  which  is  formed  in  the  first  ionic  state  4  T  has  a 

configuration  of  an  isosceles  triangle,  like  in  other  Xe2Cl 

calculations.  *  The  Xe-Cl  distance  in  the  Xe2Cl  molecule  is  3.23  A,  or 

0.16  A  larger  than  in  the  diatomic  Xe+Cl  molecule.  The  Xe-Xe  distance 

(3.25  A)  only  slightly  exceeds  the  interatomic  distance  in  the  Xe2  molecule 


(3.22  A).  The  Xe^Cl  potential  curves  for  a  fixed  Xe-Xe  distance  in  the 

Xe^Cl  molecule  (3.25  A)  are  presented  in  Fig.  4a.  The  electronic  states 

are  denoted  in  the  same  way  as  in  Ref.  30.  Because  of  the  neglect  of  the 

spin-orbit  coupling,  some  of  the  states  cannot  be  obtained  in  the  present 

calculation.  However,  in  contrast  to  Ref.  30,  the  calculation  provides  n 

states  with  wave  functions  antisymmetric  relatively  to  the  molecule  plane. 

The  number  of  states  is  3*3  *  9,  including  three  largely  neutral  states  and 

2  2  2 

six  largely  ionic  states.  Three  of  these  states  (2  4  r^,  6  are  the 

2 

states  of  n  symmetry.  The  potentials  of  the  antisymmetric  states  4  and 

2  2  2 
6  rn  almost  coinside  with  the  5  T  and  6  r  potentials  of  the  symmetrical 

states.  The  Xe^Cl  potential  curves  for  fixed  Xe-Cl  distance  of  the  Xe^Cl 

2 

molecule  (3.23  A)  are  presented  in  Fig.  4a.  The  antisymmetric  states  2  Tn> 

2  2  2 
4  rn  and  6  and  the  highest  state  8  T  are  not  shown  in  Fig. 4a.  The 

2  2  2  2  2  2 

transition  moments  for  the  pairs  of  states  1  T-4  T,  2  T-4  r  and  1  r-5  r  are 
presented  in  Fig.  4b. 

In  the  equilibrium  configuration  the  energy  of  the  Xe^Cl  molecule  is 

2 

3.09  eV.  The  dominant  transition  from  the  4  T  state  is  to  the  ground  state 
2 

1  T  with  the  transition  moment  of  2.45  D,  which  is  much  higher  than  the  ab 

28 

initio  transition  moment  of  0.96  D.  The  transition  moment  to  the  state 
2  + 

2  F  has  a  sharp  minimum  in  the  Xe^Cl  equilibrium  geometry  and  is  very  small 

2  2 

(Fig.  4b).  The  calculated  energy  of  the  4  T-l  T  transition  is  2.56  eV  (485 

nm),  which  is  close  to  the  experimental  value  of  the  center  of  the  emission 

1 1  28 
band  (2.58  eV  or  480  nm)  and  to  the  ab  initio  value  of  2.50  eV.  In  the 

DIM  calculation,"^  which  was  performed  without  taking  into  account  directly 

the  neutral-ionic  coupling,  the  transition  energy  was  found  to  be  2.74  eV. 

However,  the  difference  between  the  transition  energies  in  both  calculations 

cannot  be  explained  completely  by  the  influence  of  the  neutral- ionic 


30 

coupling,  since  in  the  DIM  calculation  the  ab  initio  diatomic  potentials 
were  used  whereas  the  present  study  applies  different,  mainly  empirical, 
diatomic  potentials. 

In  order  to  check  the  influence  of  the  coupling  term  V^.  [see  Eqs. 

(41)-(43)]  on  the  results  of  calculation,  the  Xe^Cl  molecule  was  calculated 

also  without  the  coupling  term  (V^  =  0).  By  putting  V^.  =  0  one  does  not 

neglect  the  neutral-ionic  coupling  since  the  empirical  Xe-Cl  diatomic 

potentials  include  it,  but  rather  neglects  only  the  dependence  of  the 

coupling  energy  on  the  system  structure.  The  calculation  shows  that  the 

neglect  of  the  coupling  term  V^.  decreases  the  energy  of  the  excited  ionic 

state  by  0.07  eV,  increases  the  energy  of  the  ground  state  by  0.10  eV  and, 

consequently,  decreases  the  transition  energy  by  0.17  eV,  which  becomes 

equal  to  2.38  eV.  These  energy  changes  demonstrate  that  the  influence  of 

the  coupling  term  on  the  Xe^Cl  system  does  not  change  significantly  the 

structure  of  the  Xe2Cl  electronic  states.  Such  relatively  moderate 

influence  of  the  coupling  term  on  the  Xe^Cl  can  be  explained  by  the  fact 

that  the  mixing  of  neutral  and  ionic  states  in  the  Xe^Cl  system,  as  well  as 

in  other  Xe^Cl  systems,  does  not  differ  much  from  that  of  the  XeCl  diatomic 

system.  For  example,  the  electronic  structure  of  the  Xe^Cl  molecule 

includes  3.6%  of  the  neutral  electronic  configuration,  or  1.8%  for  one  XeCl 

fragment.  The  corresponding  value  for  the  XeCl  system  with  the  same 

interatomic  distance  (3.23  A)  is  2.3%  (Fig.  2b). 

2  + 

In  addition  to  the  "main"  4  r  Xe^Cl  molecule,  there  is  also  the 

quasistable  collinear  XeCIXe  molecule,  which  has  E  state  symmetry  and  is 

100%  ionic.  The  Xe-Cl  distances  in  this  molecule  are  relatively  short,  3.06 

A  only.  The  E  state  molecule  has  0.72  eV  higher  energy  than  the  4  T  Xe2Cl 

2  +  - 

molecule  and  is  separated  from  4  T  Xe^Cl  molecule  by  a  high  barrier.  The 


energy  of  transition  to  the  ground  state  is  3.87  eV,  and  the  transition 
moment  is  3.9  D. 

2 

A  quasistable  molecule  was  found  also  in  the  ionic  4  state.  Like 
2  2 

in  the  4  T  state,  the  4  F^  ionic  molecule  has  a  configuration  of  an 
isosceles  triangle.  The  Xe-Cl  distances  are  even  shorter  in  this  molecule 

O 

(3.13  A).  However,  the  Xe-Cl-Xe  angle  is  much  larger  (85  ),  and 
consequently  the  Xe-Xe  distance  (4.23  A)  significantly  exceeds  the  Xe* 

2  4.  . 

equilibrium  distance  (3.22  A).  The  energy  of  the  4  Xe2Cl  molecule  is 

2  + 

0.8  eV  higher  than  that  of  the  4  T  Xe2Cl  molecule. 

The  ground  state-potential  energy  surface  is  found  to  have  two 

different  van  der  Waals  complexes  with  almost  the  same  dissociation  energy. 

One  of  these  complexes,  with  the  dissociation  energy  of  0.076  eV,  has  a 

symmetric  XeCIXe  collinear  configuration.  The  Xe-Cl  distances  are  equal  to 

3.24  A,  which  is  0.01  A  larger  than  in  the  XeCl  diatomic  complex.  This 

slight  increase  in  the  Xe-Cl  distance  followed  by  an  increase  of  the 

dissociation  energy  of  0.006  eV  compared  to  the  dissociation  energy  of  two 

XeCl  diatomic  complexes  can  be  explained  only  by  the  influence  of  the 

neutral- ionic  coupling.  The  second  van  der  Waals  complex  in  the  ground 

state,  with  the  dissociation  energy  of  0.075  eV,  has  an  asymmetric  geometry. 

The  Xe-Cl  distances  are  3.23  A  and  4.06  A,  and  the  Xe-Cl-Xe  angle  is  equal 

to  74.4°.  The  energies  of  the  excitation  of  the  van  der  Waals  symmetric 

2 

(collinear)  and  asymmetric  complexes  to  the  ionic  4  T  state  are  equal  to 

3.96  eV  and  3.94  eV  .respectively,  compared  to  the  excitation  energy  of 

2  2 

4.05  eV  for  the  case  of  the  XeCl  diatomic  complex.  The  1  r-*-4  T  transition 
moment  is  large  in  the  case  of  the  collinear  complex  (3.32  D)  but  small 
(only  0.71  D)  for  the  case  of  the  asymmetric  complex. 


A  van  der  Waals  complex  was  found  also  in  the  first  excited  state 

2 

(2  r).  It  has  a  symmetrical  geometry  with  the  Xe-Cl  distance  of  3.99  A  and 
the  Xe-Cl-Xe  angle  of  68.4°.  The  dissociation  energy  of  this  complex  is 
positive,  0.055  eV. 

The  equilibrium  configurations  in  the  neutral  states  were  presented 
above  as  van  der  Waals  complexes.  However,  this  term  is  not  precise  here  as 
the  attractive  forces  are  provided  also  by  the  neutral-ionic  coupling. 

B.  Xe^Cl  system 

In  order  to  check  the  possibility  of  the  attachment  of  more  Xe  atoms  to 

the  Xe*Cl  quasistable  molecule,  the  Xe^Cl  system  was  studied.  The 

formation  of  the  Xe*Cl  molecules  can  be  expected  so  far  as  the  Xe* 

molecules  were  found  experimentally.^  The  most  stable  Xe*  configuration  is 

a  symmetric  collinear  one.'’*  According  to  the  results  of  our  calculation, 

the  Xe*  arrangement  in  the  Xe^Cl  molecule  is  not  collinear.  In  the  most 

stable  configuration,  the  Cl  atom  and  two  Xe  atoms  form  an  isosceles 

triangle  with  the  Xe-Cl  distance  of  3.24  A  and  the  Xe-Xe  distance  of  3.25  A, 

2  +  - 

which  almost  coincide  with  the  interatomic  distances  in  the  4  r  Xe^Cl 
molecule.  The  third  Xe  atom  is  located  in  the  plane  perpendicular  to  the 
triangle  plane  at  a  distance  of  3.52  A  from  the  Cl  atom  and  4.41  A  from  the 
two  other  Xe  atoms.  The  energy  of  dissociation  of  the  Xe^Cl  molecule  to 
Xe*Cl  and  Xe  is  positive,  0.113  eV.  It  follows  that  the  Xe^Cl  can  be 
formed,  although  the  Xe^Cl  -Xe  binding  is  weak  and  is  determined  mainly  by 
the  polarization  interaction.  The  third  Xe  atom  bears  only  0.5%  of  the 
positive  charge.  The  transition  energy  to  the  ground  state  is  2.48  eV. 

This  value  is  not  far  from  the  Xe^Cl  transition  energy  (2.56  eV),  which 
makes  it  difficult  to  distinguish  the  Xe*Cl  and  Xe^Cl  molecules  by  their 
emission. 
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A  less  stable  Xe^Cl  molecule  is  formed  in  the  symmetric  configuration 
with  the  Xe-Cl  distance  of  3.33  A  and  Xe-Xe  distance  of  3.58  A.  The 
dissociation  energy  in  this  configuration  is  0.08  eV,  and  the  transition 
energy  is  2.62eV. 

C .  Xe^Cl  system 

According  to  the  results  of  the  calculation,  the  spherical  symmetric 
Xe^Cl  system  is  unstable.  The  stable  configuration  is  formed  by  a  Xe^Cl 
molecule  with  slightly  extended  interatomic  distances  (3.29  A  for  Xe-Cl  and 
the  same  for  Xe-Xe)  and  10  weakly  bounded  Xe  atoms.  These  10  Xe  atoms  bea 
together  only  4.3 %  of  the  positive  charge.  The  energy  of  the  dissociation 
of  this  system  to  Xe^Cl  and  10  Xe  atoms  is  0.88  eV.  The  energy  of 
transition  to  the  ground  state  is  2.49  eV.  Considering  the  dissociation 
energies  of  the  Xe^Cl  and  Xe^Cl  systems,  we  have  to  note  that  they  may  be 
overestimated  because  of  the  neglect  of  three-  and  four-center 
interactions."^ 

VI .  Conclusion 

1.  The  diatomic-in-ionic-systems  (DIIS)  method  is  disigned  for  a 
semiempirical  calculation  of  neutral  polyatomic  systems  R^X,  which 
consist  of  atoms  R  of  S  symmetry  and  an  atom  X  of  P  symmetry  and  which  can 
form  excited  ionic  R*X  states  with  positive  ions  R+  of  P  symmetry  and  a 
negative  ion  X  of  S  symmetry.  Rare  gas(R)-halogen(X)  systems  satisfy  this 
definition  and  have  been  subject  of  consideration  in  the  present  paper.  The 
DIIS  method  can  be  also  used,  with  some  modifications  in  the  matrix  element 
expressions,  for  calculations  of  other  systems,  such  as  hydrogen-containing 

rare  gas-halides  R  H  X,  alkali  metal-halides  and  rare  gas  cluster  ions  R+. 

8  n  m  n 

2.  The  DIIS  method  allows  us  to  calculate  the  electronic  structure  of 


polyatomic  systems  by  using  diatomic  potentials  and  transition  moments  as 


initial  data.  In  order  to  improve  the  results,  it  is  possible  also  to  use 
some  known  polyatomic  values  as  semiempirical  parameters.  However  the 
present  calculation  of  polyatomic  systems  is  based  on  diatomic  data  only. 

The  calculation  is  limited  to  states  without  spin-orbit  coupling.  The 
influence  of  this  coupling  on  the  system  energy  is  taken  into  account 
approximately  and  indirectly  by  using  empirical  diatomic  potentials. 

3.  The  DIIS  method  has  been  applied  to  the  Xe^Cl  systems,  Xe^Cl  in 

particular.  The  ionic  Xe^Cl  in  its  most  stable  state  (4  T)  was  found  to 

have  an  equilibrium  geometry  of  an  isosceles  triangle,  like  in  other 
28  30 

calculations.  ’  In  the  equilibrium  configuration  the  dominant  transition 

is  to  the  ground  state,  although  in  other  configurations  the  radiative 

transition  to  the  first  excited  state  becomes  possible  as  well.  In  addition 
2  + 

to  the  4  T  Xe^Cl  molecule,  there  are  other  quasistable  ionic  molecules  with 
higher  energy,  in  particular,  the  collinear  Xe^^+Cl  Xe*^+  molecule.  Any 
empirical  evidence  of  the  emission  of  these  ionic  molecules  are  unknown  to 
us.  In  the  ground  state,  two  van  der  Waals  complexes  were  found  with  the 
sane  dissociation  energy  but  different  geometries  and  different  transition 
moments  to  the  ionic  states. 

4.  According  to  the  results  of  the  calculation,  the  ionic  Xe^Cl 
molecule  attracts  more  Xe  atoms,  but  without  sharing  with  them  much  of  the 
Xej  positive  charge.  The  ionic  molecules  with  more  than  two  Xe  atoms, 

Xe^Cl  in  particular,  can  be  formed  in  a  dense  gas.  It  is  difficult  to 
detect  these  molecules  as  their  emission  wavelength  is  similar  to  that  of 
the  Xe^Cl  molecule. 

5.  The  symmetric  Xe^Cl  clusters  were  found  to  be  unstable,  so  that 

13  20 

their  formation  in  rare  gas  solids  '  is  unlikely. 
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Table  1 


Parameters  of  diatomic  terms.  For  the  designation  of  the  parameters  see 
Eqs.  (44)  and  (45).  Rg  and  De  are  the  equilibrium  distance  and  the 
dissociation  energy,  respectively. 


uo 

a 

A 

C6 

C4 

C1 

R 

e 

D 

e 

diatomic 

fragment 

term 

eV 

A*1 

eV 

eV  A6 

eV  A* 

eV  A 

A 

eV 

Xe-Xe 

I 

0 

3.17 

24170 

320 

0 

0 

4.47 

0.023 

Xe+-Xe 

UI 

0 

2.95 

14520 

320 

29.14 

0 

3.22* 

1.08* 

(Xe-Xe)+ 

VI 

- 

1.45 

-172 

- 

- 

Xe+*Xe 

un 

0 

3.00 

18625 

320 

29.14 

0 

4.00* 

0.16* 

(Xe-Xe)+ 

vn 

0 

1.90 

-168 

- 

- 

XeCl* 

z 

0 

3.85 

31670 

60 

29.14 

0 

3.30 

0.195 

8.0** *** 

4.108 

Xe-Cl 

11 

0 

3.10 

1167 

99.2 

0 

0 

3.23 

0.035 

Xe+-Cl" 

22 

8.52 

2.512 

2033 

0 

66 

14.4 

3.01 

hSs 3 

Xe-Cl 

in 

0 

3.90 

86000 

99.2 

0 

0 

4.10 

0.015 

*** 

310 

4- 

‘f.jc 

Xe-Cl 

2n 

8.52 

2.847 

5520 

0 

66 

14.4 

3.074 

3r97 

Xe-Cl 

vi 

- 

1.20 

-30 

- 

- 

- 

- 

- 

Xe-Cl 

vn 

- 

1.40 

-18 

- 

- 

- 

- 

- 

*  Re  and  D@  for  Xe^  potential  curves  connected  with  U  and  V  by  Eqs. 
( 35 ) - ( 38  ) . 

**  Parameters  of  the  second  exponent  of  the  XeCl  potential. 

o 

***  Coefficient  Cg  of  the  Cg/R  term  of  the  III  potential. 
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Figure 
Fig.  1 
Fig.  2 


Fig.  3 


Fig.  4 


Captions 


XeCl  potential  curves  (Table  1). 

The  effects  of  coupling  between  ionic  and  neutral  electronic 
configurations  in  the  XeCl  system.  The  vertical  dashes  stand  for 
the  211  (C)  equilibrium  distance. 

a)  Stabilization  energy  of  the  coupling  (right  ordinate)  and  the 
contribution  of  ionic  configuration  (1 ««  ordinate)  in  the  ground 
II  (X)  state. 

b)  1Z-2Z  (X-C)  and  111-211  (A-C)  transition  moments. 

Xe^Cl  potential  energy  surfaces  for  a  geometry  (shown  in  the 

figure)  with' a  fixed  Xe-Xe  distance  equal  to  that  of  the  Xe^Cl 
2 

(4  T  state)  molecule.  The  vertical  dashes  stand  for  the 
equilibrium  Xe^Cl  configuration. 

Xe^Cl  system  in  a  geometry  (shown  in  the  figure)  with  a  fixed 
Xe-Cl  distance  equal  to  that  of  the  Xe2Cl  (4  T  state)  molecule. 
The  vertical  dashes  stand  for  the  equilibrium  Xe^Cl 
configuration. 

2 

a)  potential  energy  surfaces  for  five  r  states. 

b)  l2r-42r,  l2r-52r  and  22r-42r  transition  moments. 
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